A Non-Extensive Model for Quark Matter Produced in Heavy Ion Collisions 



Tamas S. Biro and Gabor Purcsel 
MTA KFKI Research Institute for Particle and Nuclear Physics, H-1525 Budapest, P.O.B. 49, Hungary 

We describe quark matter in the framework of non-extensive thermodynamics. We point out 
that particle spectra with power-law tail lead to an increased energy and entropy per particle, and 
therefore even a massless plasma may fit both the relatively high value of E/N — 1 GeV and the 
spectral slope of T = 175 MeV observed in RHIC experiments. 
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^ ■ I. INTRODUCTION 

^ Cy The probably most important motivation for relativistic heavy-ion experiments at large scale accelerator facilities 

^ (RHIC, SPS) is to study quark matter in a unique environment, possibly in form of a quark-gluon plasma (QGP) 
produced for a short while before hadronization While the original idea of quark-gluon plasma sticks to an 

' equilibrium picture both in its early, bag-model inspired versions 0| and in the theoretically most advanced lattice 

^ ' gauge models 's', , phenomenological studies of the experimental spectra offer an increasing amount of evidence that 

^ ' this equilibrium is not completely achieved IH.I3- 
00 

In particular the equilibrium thermodynamics relies on (in the Boltzmann limit) exponential particle spectra, while 
experiments definitely show a power-law tail at high transverse momenta. This means that some assumption leading 

, to the Boltzmann-Gibbs factor is not fulfilled. We have therefore to look for a theoretical description which is more 

' general, which releases one or more original assumptions. Such a candidate is the non-extensive thermodynamics, 

[ promoted by Tsallis and others, which is based on an altered definition of the Boltzmann entropy P. ll0lllllir2, 13, 14|- 

• More recently it has been shown that the power-law distribution also can be described by a logarithmic dispersion 

[ relation for the quasi-particle energy while keeping the original, extensive entropy definition [l^l . Il(i| . In this case 

^p-i. the energy is non-extensive. The relevance of Tsallis statistics to high-energy e+e~ and heavy-ion collisions has been 

O i' studied in several papers |17| . In this paper we show that such a version of non-extensive thermodynamics is able to 

, describe a power-law distribution of quasi-particle energies and explain an apparent spectral temperature, T — 175 

' MeV at low energy. Furthermore, the very power in the tail of this distribution can be connected to the energy per 

^ , particle, E/N = 1 GeV, which remained unexplained in traditional thermal models. 



X 



There were several attempts to describe dynamical mechanisms which would simulate (tot ally or just in part) a Gibbs 
distribution or a Tsallis distribution with a power-law tail asymptotics [iM [TqI l20l l2ll |22 | . Tsallis advocates the 
non-extensive thermodynamics based on his entropy definition for a general use for systems in non-ergodic states 
showing exponentially not suppressed distribution tails (2^ (2j| . Due to this an interesting question arises for the 
heavy ion physics: can the quark-gluon plasma (QGP) be described in terms of the non-extensive thermodynamics? 
Would it give a quantitatively better description of experimental observation than the old bag-model picture while 
being almost as simple? In this paper we shall answer these questions positively by showing that a simple QGP 
with power-law energy distribution instead of the Gibbs one and a bag constant is able to reach the hadronically 
fitted temperature of Tq w 175 MeV and the 1 GeV energy per particle at the same time. We extend this study 
for finite chemical potential using relativistic Fermi and Bose distributions derived according to the rules of non- 
extensive thermodynamics. After briefly discussing the recently considered directions in this research of extending 
thermodynamics, we choose the interpretation of Q. Wang [ia,[l9> who has revealed that an anomalous (logarithmic) 
quasi-particle dispersion relation may lie in the background of the canonical Tsallis distribution with the famous 
power-law tail. We shall study the equation of state (eos) of this QGP and eventually obtain a constant energy per 
particle curve on the T — /i plane. 
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II. QGP IN NON-EXTENSIVE THERMODYNAMICS 



The starting point of our description of a quark-gluon matter at the instant of hadronization is based on a power-law 
energy distribution of partons, 

w^{E) = (l + Elbr\ (1) 

containing two parameters; an energy scale h in the order of 1 to a few GeV and the power c, which is expected to 
be rather high, 4-10, extracted from fits to minijet distribution in a pQCD based theoretical parton dynamics model 
[25I l2^. The ZEUS collaboration extracted the value c — 5.8 ± 0.5 in e+e~ experiments j23|. In the double limit, 
c ^ 00, & ^ 00 with a fixed ratio T = b/c this distribution approaches the familiar canonical Gibbs-Boltzmann 
formula, 

«;oo(-B) =cxp(-£;/T). (2) 

As it has been indicated in Ref.|0| parton recombination, which is most likely to happen at low relative energies and 
therefore combines a parton cluster or pre-hadron with energy E from two partons with energy E/2, or in general 
from n partons each with an energy of E/n^ drives the power-law distribution towards the exponential limit: 

Wn{E) = {\ + El{nb))-'^^, (3) 

by increasing both the effective c and h to nc and nh while keeping their ratio T — nb/nc = b/c constant. The 
power-law distribution we are using here is related to a canonical Tsallis distribution by setting the temperature 
T = b/c and the Tsallis index tog = 1 + 1/c when considering w'^-weighted averages of energy and particle number 
as thermodynamical observables. Funny enough, using the "normal", w- weighted averages as macroscopic variables 
one arrives at q = 1 — 1/c. 

An anomalous thermodynamics based on an altered definition of the entropy, 

i 

has been worked out by Tsallis and several other authors 0, 0, 0, 0, 0| . It resembles the Legendre transforma- 
tion structure and herewith the canonical and grand canonical versions of phase space distributions behind observable 
macroscopic expectation values. The distribution Q can be viewed as the canonical Tsallis distribution 

1 / (q~ 1)E,\^ 
-. = ^(l + ^S^) (5) 

with the canonical partition function. 

The interpretation of the temperature T and the Tsallis-index q (related to the power-low cutoff b and power c) 
is still subject to the fact, however, that which version of macroscopic observables are considered in the effective, 
"non-extensive" thermodynamics. By now the choice 

E = J2^lEjJ2^l (7) 

i i 

and the canonical variational principle, 

S-PE = max. (8) 

seem to establish themselves. (However, a recent article jl^ righteously criticizes this choice. According to the 
interpretation given there the "normal" macroscopic averages are meaningful as well.) In the grand-canonical version 
E is replaced hy X = E — iiN, and Ei by Xi = Ei — fiQi. The expression to be minimized becomes S — (3X. 
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In the following we review the most important details of these calculations. Considering the derivatives with respect 

to Wi, 

—S = cil-qwJ), 

d w.' 

Ti — ^ = 
dwi 2^w 

one arrives at the constraint 



X = ,^iXt-X), (9) 



This gives rise to a grand-canonical distribution wt proportional to a (— c)-th power: 



This can be casted into the Tsallis form, 



(12) 

normalized by the partition sum, 

The temperature, T is, however, no more related to the Lagrange multiplier /? = 1/6 as straight as usual. It is 

T = e"^w]-X/c. (14) 

Since the original probability is normalized, 

^Wi = l, (15) 

i 

in general J^i differs from one. 

In order to make the discussion more transparent we use the following abbreviations for the high-power-law distribution 
function and its inverse: 

exPc(a;) =(l-x/c)-^ — >cxp(a;), 
Inc(x) =c{l-x-^/'') — >lii{x). (16) 

The limits are valid for c — > oo by fixed T = b/c. These functions are inverse to each other, expg(lnc(x)) = x, and 
reflect the non-extensive property, 

luc {xy) = luc {x) + luc (y ) - i luc (c) luc {y). (1 7) 

The derivatives of these functions follow rules similar to the familiar ones: 

-^exp^(a;) = exp«(a;), 
^lne(x) = 1/x'^. (18) 

The Tsallis entropy with this notation reads as 

5 = ^u;ilne(lM), (19) 
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and in case of equiprobability, Wi = 1/N for N states, one arrives at 

S' = lnc(iV), and N ^ exp^{S). (20) 
The grand-canonical distribution and partition sum becomes 

= I exp,(-X,/r), Z^Y. exPc(-^./r). (21) 

i 

The anomalous thermodynamics reflects famihar looking relations, with some exceptional points. The heat is expressed 
by 

esj^wl = T\ncZ + E- fiN, (22) 

i 

with the grand canonical partition sum 

Z = Y,e^Pci-iE^~f^Q^)/T), (23) 

i 

where we use the q- weighted particle number, 

(24) 

i i 

and the chemical potential, /i, the corresponding Lagrange multiplier. Here Qi are the conserved charges of the state 
i. The pressure, assuming homogeneity in a sufficiently large volume, is given by 

p = ^ In, Z. (25) 

Resolving expression H14|l for O one arrives at the following grand canonical potential 

- n = pV = Thic Z = TS - {E - fiN){l - S/c). (26) 

The derivatives are given by 

^^-~P^ = ^n^-N{l~S/c). (27) 

For the energy per particle the correction factor (1 — S/c), by chance, does not play a direct role: 

J, IN Kdp/d^i)+T{dp/dT)-p 



In case of an ideal (but still power-law, not Gibbs-distributed) gas one assumes that 
with 



I TT^Z^TWZ (29) 



Z = exp,(— nx), (30) 

n 

and 

Xk = Ek- ^lQk■ (31) 

(Here Qfe = Q, the conserved charge is usually independent of the momentum of the considered particle.) For fermions 
the summation runs over n = and n = 1 only, while for bosons for all positive integers n = 0, 1, 2, . . .. In this paper 
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we shall deal with massless relativistic particles (quarks and gluons), so practically Ek = \k\. The derivative relations 
known for the ordinary ideal gas are almost fulfilled, 

n=— , e = /in + es-p (32) 

for the particle-, entropy- and energy-density, respectively. 

Constructing this way the Tsallis-Bose and Tsallis-Fermi distributions (alternative ways can be found in 28]), one 
may use the following integral representation 4^1 of the power-law function: 

exp,(-:E) = — - / dtt^-ie-*(i+"/^) (33) 
r(c) Jo 

where r(c) is Euler's Gamma function giving (c— 1)! for integer c > 1. This representation relates the grand canonical 
partition function Z of anomalous thermodynamics to the usual one-particle forms of Fermi and Bose distributions, 
respectively: 

where C(y) = 1-1- exp(— y) for the Fermi, and C(y) = 1/(1 ~ Gxp(— ?/)) for the Bose distribution. This form is 
quite suggestive: the Tsallis thermodynamics is an average of "normal" systems with an Euler-Gamma distributed 
parameter, factorizing the energy. Applying it for x — E/T we substitute y — tE/b, the scaled energy into the 
"normal" partition function (.. This sheds some light on interpretations of Tsallis-distributed particle spectra in terms 
of a temperature fluctuating according to the Gamma distribution |2lj | or energy non-conservation due to higher 
than two-body collisions, again with Gamma distributed energy imbalance |29l l30j|. These correspondences are, of 
course, no explanations for the power-law distribution, but rather equivalent re-formulations, assuming the Gamma 
distribution for the corresponding (Legendre associated) intensive variable. 

Our main goal in this paper is to investigate the original idea of quark-gluon plasma in terms of the anomalous, non- 
extensive thermodynamics. Therefore we do not wish to speculate here about the microscopic dynamical mechanisms, 
which lead to the power-law distribution, closely approximating exponential spectra in the intermediate transverse 
momentum range. We assume now, that an anomalous QGP of massless and Tsallis distributed particles have existed 
before hadronization and ask the question whether this state is realistic, whether it is nearer or farer from the state of 
hadronic matter observed in RHIC (and partially in SPS) experiments with many of its properties surprisingly well 
fitted by traditional equilibrium thermodynamical concepts underlying the hadronic thermal model |3lL 133 . 133 . l34| . 
For this purpose we select the energy per particle E/N as a decisive observable; it can be quite directly measured 
on the one hand and quite easily calculated in thermodynamical framework on the other hand. Considering a bag 
constant also in the case of Tsallis QGP, we shall modify E/N accordingly. Another interesting quantity is the entropy 
per particle, S/N. It can be related to other thermodynamical quantities as 

l = l = '-±P-!^. (35) 
N n riT T 



III. ENERGY PER PARTICLE 



For calculating E/N it is helpful to introduce notations for the (1 — S/c) uncorrected derivatives of the pressure, 
Co = Tdp/dT -t- ^dp/d^ — p, and Uq — dp/dfi. Since the derivatives of the pressure with respect to T and fj, can 
be regarded as derivatives of Inc Z{x) with respect to x under the phase space integral by noting that d/dx acts as 
—Td/dEk or Td/dji respectively, we get 

f (Pk Ek - ^J,, 
"° = y (2^^=(^^)' 

- /(g,i^./c(^), (36) 
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with the one-particle distribution function 



Now we have simply E /N — eo/uQ. 



fcix)^-^lncZ{x). (37) 



Before dealing with Fermi and Bose statistics in the anomalous thermodynamics let us first review the dilute matter 
limit, the Tsallis-Boltzmann approximation. This relies on the assumption that exp^(— x) ^ 1. Now Z = \ to leading 
order and there is no difference between Fermi and Bose distributions, 

/,(x)«--^lne(H-exp,(-x))«--^exp,(-a:) = + '\ (38) 
ax ax \ c/ 

A (slight) difference can be observed between the inverse slope of logarithmic one-particle spectra and the temperature, 

(c + l)r3iopc = cT = b. (39) 

We carry out our first estimate at fi — and finite temperature in the Boltzmann approximation. Using fc{E) = 
(1 + E/b)^'^-^^ in the corresponding integrals we get 

d 2b^ 

no 



27r2 c(c-l)(c-2)' 
_ _d 6b^ 

" 27r2c(c-l)(c-2)(c-3)' ^ ' 

with a common (color, spin, etc.) degeneracy factor d. In the low energy region we obtain a simple relation between 
the energy per particle and the inverse slope of the one-particle spectra, 

E/N - eo/no - ^3T = ^lili^T^iope- (41) 
c — 3 (c — 3) 

This differs from the familiar relation, E/N = 3T, known for the massless Boltzmann gas. It is interesting to note, that 
- as it has been calculated in Ref. |35ll3B | - by using the normal weighted definition for particle number, N' — J2i WiQi, 
one gets E/N' ~ 3r'495. We do not think, however, that using different definitions for the macroscopic energy and 
particle number would make a sense in comparison with experimental data: either Wi or wf is the distribution of 
physical states. Furthermore T is not exactly the "experimental temperature" conjectured from the inverse spectral 
slope. 

An estimate of the inverse slope gives rise to Tgiopc ~ 175 MeV for RHIC experiments in the /i ~ region. This 
is also the value used by conventional thermal models for the temperature. Aiming at E/N = 1 GeV now, the 
formula H41I) leads to a pow er of c-l- 1 « 8.42 in the one-particle distribution fc{E), a value comparable with minijet 
PT-distribution fits [25l l26l|. Neither the traditional Boltzmann formula, E/N = 3T, nor the simple bag model of 
QGP, where E/N = 4T, is able to describe both this inverse slope and energy per particle at the same time. 

In order to display a more throughout comparison to the thermal model points in the T — ^ plane fitted to different 
accelerator experiments (47j , we carried out numerical integrations with Tsallis-Fermi and Tsallis-Bose distributed, 
massless quark-gluon plasma at finite temperature and chemical potential. In this case, denoting a general integral 
of fc{x) weighted with a power x" by 

/oo 
dxx^Ux), (42) 
-a 

and the symmetric and anti-symmetric combinations, Snia) — ^n(ti) + ^ni^ci)^ ^n(^) ln{0') — In{—ci)-, we arrive at 

27r2no = T^S2{^i/T) + 2T^^lA^{^x/T)+T^l^So{^l/T), 

2^2 eo = T^Ssip/T) + 3TV^2(M/r) + Sli^^/T) + T^i^Aoi^i/T), 

(43) 
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for each fermion. For the boson part we assume /i = 0. The energy per particle characteristic for the power-law 
distributed mixture then becomes 

^Ij^ ^ dseoBiT) + dpeoFiT, fi) + dpepAFiT, fi) ^^^^ 
dsnoBiT) + dpnopiT, ^i) + dpnoApiT, n) ' 

Before presenting our numerical results there is one more aspect of this version of the anomalous thermodynamics 
to be discussed. In the traditional, Gibbs thermodynamics the anti-fermion distribution is related to the fermion 
distribution simply by a change in the sign of the baryochemical potential, fi. This formula exactly coincides with 
the density of fermion holes: 

l-foo{i-E-fi)/T) = f^{{E + fi)/T) (45) 

due to foo{x) + foo{—x) = 1 for the Fermi distribution foo(x) — 1/(1 + exp(a;)). This relation is not fulfilled for a 
finite c index Tsallis-Fermi distribution: 

/e(a;) = (l + exp_,(x))-«. (46) 

In this case 

Mi-E - ii)/T) + f,{{E + ii)/T) - /e(-x) + f,{x) ^ 1, (47) 

instead 

fl''^{x)+fj,\-x) = l (48) 

is fulfilled. Here the anti-charge choice with —jj. and the particle hole choice by replacing —E energy are not equivalent. 
For the estimate of the energy per particle of a Tsallis QGP at finite temperature T and baryochemical potential /i 
in Ref.Q we utihzed the particle-hole correspondence: 

7,(i?;M,T) = l-/,(-i?;^,T). (49) 

Now we turn to another interpretation of the non-extensive thermodynamics, similar to that suggested by Q. Wang 
|l5lll6ll3^ . which handles the fermions and antifermions in a more satisfactory manner. 



IV. QUASIPARTICLES AND ANOMALOUS THERMODYNAMICS 



In the preceding discussion we have already indicated some disadvantages of the " main stream" Tsallis formulation 
of non-extensive thermodynamics. These make it to be not the best choice for our purpose, i.e. for using it in the 
description of quark matter. The dependence of some thermodynamical relations indirectly on the system size (i.e. 
dependence of densities on the total entropy, etc.), the slight difference between temperature parameter and one- 
particle distribution inverse slope, and finally the strange correspondence between particles and holes in the Fermi 
distribution leave us with dissatisfaction. Also the very basic point, the non-extensivity of the Tsallis-entropy remains 
a challenge for the physical interpretation (and acceptance) of this approach. 

Fortunately it has been shown meanwhile that a mathematically well-defined monotonous function of the Tsallis 
entropy describes an extensive entropy measure jlSj, ,3SJ . The Tsallis-entropy satisfies the general rule of pseudo- 
additivity, 

/(^12) = /(5l) + f{S2) + A/(5i)/(52), (50) 

required by the existence of thermal equilibrium in non-extensive composite systems. The Tsallis-entropy realizes the 
case f{S) — S and \ — \/ c = q — 1. The additive, i.e. extensive measure of entropy then becomes 

^ = i ln(l + A/(^)). (51) 

Based on the maximalization of S of the Tsallis-entropy instead of S itself, a formally traditional thermodynamics and 
canonical, or grand canonical distribution emerges The price to pay for that is, that the energy (and incidentally 
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the particle number and the volume) is no more an extensive quantity. It is a clear indication that the Tsallis entropy 
describes systems whose parts still interact with long range forces. The Tsallis canonical distribution and partition 
sum can be transformed into an ordinary canonical partition sum using a " deformed Hamiltonian" . By defining 

H - ^) In (1 + (1 - q)m - ^iN)) (52) 



one arrives at 



Z = Tr (f - (g - 1)(3{H - ^iV))^/('-^^ = Tre"''^. (53) 



As a consequence the Fermi and Bose distributions keep their traditional forms as functions of the eigenvalues of the 
deformed Hamiltonian, H. The trace over the states of the many particle system leads to a product over the phase 
space labels of independent quasi-particles with energy ujk, and the one-particle distribution becomes 

where E^. is the eigenvalue of the original Hamiltonian, i.e. the energy of an interacting quark or gluon with momentum 
k. For particle and nuclear physics purposes this use of a quasi-particle energy is the more enlightening approach. 

For the investigation of a long range interacting quark matter we still have to modify the above picture a little. The 
power-law tail with a general, fractional power restricts the base of this power-law to be positive. In a picture of a 
QGP rooting in QCD on the other hand gluons are exactly and quarks are practically massless, therefore Ek = \k\ 
can be as small as zero. For any finite, positive chemical potential /x > the above distribution may become complex, 
unless fi < cT = b. This restricts the use of quark chemical potentials too seriously. We explore therefore here another 
approach. 

Repeating the above line of derivation given by Q.Wang we start with the anomalous dispersion relation for the QGP 
quasi-particles: 

iJk = b\nil + Ek/b) (55) 



with Ek — yk^ + rr? being the dispersion relation of a free relativistic particle with mass m. FigQ] shows this for 
different masses. The right side of the figure plots the group velocity, Vg = duj/dk, which overall remains below one. 
Such a dispersion relation does not violate causality. Still it is anomalous; it shows an acoustical branch for the 

massless case, and a line crossing from the optical branch at low momenta to an acoustical behavior at high momenta 
for finite mass. Such a quasi-particle dispersion relation has never been considered yet in QCD, as far as we know 
[mii^lisf . The plasmons all live on the optical branch, above the light cone. Below the light cone a series of Landau 
poles can be located, they cannot belong to traditional " particles" . It is also hard to find a (partial) summation of 
diagrams leading to a logarithmic dispersion relation from QCD. But the proof of impossibility is also unknown 50]. 

The investigation of the non-rclativistic and extreme relativistic limit of the formula H55|l may help to approach to a 
physical interpretation. Based on the second derivative, 

d'^uJk ^ mVtjfc 4- (m^ - fc^)/6 

dk^ " (LOk+Uol/b)^ ^ ' 

we get in the non-relativistic range <C fc ^ m (which nevertheless cannot be applied to massless QGP) 

^fc = 61og(l + -) + ——— + ... (57) 
+ m Zm 

This features an attractive mean field effect for the heavy quarks in QGP due to luq < m, and an increased effective 
mass meff = m{l + m/b). For relativistic (including m — 0) particles the quasi-particle mass cannot be interpreted 
other than the inverse of the second derivative of the dispersion relation at = 0, which is zero as defined above. At 
any finite k the second derivative is negative, not corresponding to any picture of a traditional point particle. 

It is also interesting to interpret this dispersion relation in the framework of plasma physics. In plasmas a dielectric 
coefl[icient, e(w, fc), may alter the free dispersion relation. The pole of the propagator in medium is located at 

D-^ =uj^ - k^/e = 0. (58) 



1 2 3 4 5 

k/b 



1 2 3 4 5 

k/b 



FIG. 1: Anomalous dispersion relations and group velocities for Tsallis quark matter. 



From this we get 



^ ujk ln(l + fc/6) ^ ' 



which is greater or equal to one. This is a further anomalous property, related to a negative self-energy (due to 
Z?~^ = up' — — li, li = fc^(l/e — 1) ). This property also occurs in non-abelian gauge theories, where gluon loop 
contributions render the gluon self energy negative. This effect is alike asymptotic freedom, a vanishing charge at 
high momenta. Classically a = ao/e, according to ea. H59(l this quantity tends to zero at fc ^ oo. The very functional 
form of this running coupling, however, differs from the usual 1/log one, known from perturbative QCD. 

The one-particle distribution functions (Fermi and Bose) in the Boltzmann approximation become simple power-laws: 

/ ^ \-b/T 

f{k) = cxp(-(c^fe - ^i)/T) = cxp(Ai/r) f 1 + -fj (60) 

with a power c — b/T. Now the chemical potential /i does not enter under a logarithm, which makes it possible to 
obtain real results at any value of it. The Fermi and Bose distribution becomes 

^^^^ " exp(wfc-M)/r±l " '^^^Pi^jTW+EiJl)f±\' ^^^^ 

In the c — > oo limit the equilibrium thermodynamics of an ideal gas is re-established. Otherwise a power-law tail 
occurs for momenta around and higher than b. Fig|21 shows these distributions for the value b = 1.0325 GeV, favored 
by the RHIC experiment slope T — 175 MeV and an energy per particle of E/N — 1 GeV. 
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FIG. 2: One-particle distributions in Tsallis-QGP. Quarks or gluons correspond to c = 5.4496, mesons to 2c and baryons to 3c 
according to the recombination idea. An exponential fit to fc = is also indicated. 

This anomalous behavior can be traced back to one feature: the familiar exponential equilibrium distribution has 
been modified at high momenta. This contradicts to the traditional result from perturbative QCD, so it cannot stem 
from an interaction correction to the static, equilibrium state at finite temperature. We rather think it may be related 
to a characteristically short time scale preventing equilibrium. The short-time behavior is somehow non-particle like, 
probably related to high off-shell-ness of emerging hadrons. 

Recovering the power-law in the Tsallis-Boltzmann approximation, the important thermodynamical quantities coincide 
with our previous results [ref RHIC]. An other important limit, the T — Fermi distribution also can be calculated 
analytically for massless quark matter using the logarithmic dispersion relation (|55() . The contributions to pressure, 
energy density and particle density follow the traditional relation, e = Ts + fin — p. Considering a bag constant, 
e — cq + B and p = pq — B still satisfy this relation. The stability edge line, where p = 0, can be realized by po — B, 
and therefore e = cq +po- On this edge of the QGP stabihty against clustering, - possibly a state at the hadronization 
process, - the energy per particle and the entropy per particle are linearly related for any ideal quasi-particle system. 
In particular at T = 0, E/N = /i has to be fulfilled. As a consequence a hadronic thermal model with E/N — 1 
GeV would have to end at /is = 1 GeV. This estimate does not leave much room for /i = /is/S massless quarks as 
immediate precursors of the hadrons at the same energy per particle at T = 0. 

Direct integration under the Fermi energy /i gives the following results at zero temperature: 





(62) 



n 



with X = kp/b. The Fermi momentum can be obtained by inverting the logarithmic dispersion relation, kp = 
^(gM/fc _ \Y Considering baryons made of three quarks, finally we set /Lts = 3/x. In the b ^ oo limit, realizing that the 
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FIG. 3: Energy per particle in Tsallis-QGP without a bag constant. 



pressure subtracts from the first three terms in the Taylor expansion of ln(l + x) exactly ln(l + x), it is easy to see 
that po = (247r^). In this small x (kp ^ b, dilute, cold quark matter) Hmit eo/n ~ 6(3x/4 + 9x^/10 + .t-^/G + . . .), 
so the pure energy per particle is eo/n = 3/z/4. 



V. NUMERICAL RESULTS 



Finally we present results of numerical integrations at arbitrary T and with logarithmic dispersion relation for 
massless quark matter. First, FigOlplots the direct energy per particle eo/n as a function of temperature for different 
chemical potentials, utilizing the h = 1.0325 GeV value. The requirement of having E /N = 1 GeV at /i = in this 
case leads to T = 195 MeV. Alternatively T = 175 MeV is achieved by the choice of 6 = 0.8275 GeV. 

The direct eo/n = 1 GeV lines on the T — ^ plane are shown in Fig0|The constant power, c = 4.72732 line (solid line) 
approaches the data of the hadronic thermal model (boxes) most closely, however only with the /i = assumption. 
This is not readily a quark matter, where ii = ^b/3 would be. Massless baryon clusters as prehadrons satisfy rather 
this assumption, but this idea - however exotic on its own right ~ does not have yet any support from microscopic 
models. The constant dispersion relation with b = 0.82728 GeV plots another line. This deviates even more at T = 
from the hadronic boxes, contrary to the fit to the same point at fi — 0. (Both the c and the b value were obtained 
from E/N — 1 GeV at = and T = 0.175 GeV.) Finally the traditional, massless ideal gas corresponding to the 
exponential spectra (c = oo) are vastly off from the experimental data. 

The entropy per particle is a sensitive observable, because cannot decrease in spontaneous hadronization processes 
unless there are more hadrons than quarks and gluons before. It is plotted in the massless ideal Tsallis-gas scenario 
along the E/N = 1 GeV curve as a function of the temperature. The solid line corresponds to a constant c = 4.72732. 
It is worth to note that it ends at a finite value at T = 0, contradicting to the third law of thermodynamics. Probably 
the c =constant assumption is not realistic. (Of course this deviation is proportional to 1/c as it is easy to see from a 
low-T expansion of the Fermi distribution.) The constant b = 0.82728 GeV curve (dashed-dotted) and the traditional 
thermodynamics (dashed) curve both start from zero at T = 0. However, only the finite b case can fit the high value 
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0.35 




H (GeV) 

FIG. 4: The 1 GeV energy per particle line without bag constant in the T ~ fi plane. The constant power, c = 4.72732, line is 
solid, the constant energy parameter, b = 0.82728 GeV line is dashed-dotted, the traditional thermal line (c = oo) is dotted. 
Boxes show corresponding values from the hadronic thermal model. 

of S/N w 6.7 at T — 0.175 GeV dictated by the E/N = 1 GeV for massless bosons and fermions. 

Taking into account a possible bag constant in Tsallis-QGP with logarithmic " quasi-particle" energy, the fitted b 
value becomes 1.0325 GeV. In FiglHIthe (kinetic + bag) energy per particle, (eo +pa)/n at the p — stability limit 
is plotted as function of the temperature, T (left part), and as a function of the chemical potential, /jb = Sfi (right 
part). 

Finally we present the p = stability line for a constant b = 1.0325 GeV dispersion relation, while adjusting the bag 
constant T and n dependent so, that the total energy per particle, E/N = (eo +Po)/n remains constant, 1 GeV. 
We find that of all Tsallis-distributed massless QGP models this prescription may come the closest to the hadronic 
thermal model data. Still the low value of the hadronic fits to exponential spectra cannot be achieved, pointing 
out that massless QGP at lower-energy heavy-ion experiments is probably not the state of matter at hadronization, 
not even with power-law high-momentum tail. Instead massive quarks and gluons have to be considered, which would 
allow to reach high E/N and S/N values with high-pt powers closer to the experimental values. At RHIC energies, 
however, this picture is promisingly consistent. 



VI. CONCLUSION 



In conclusion we have worked out the thermal properties of a plasma of massless quarks and gluons utilizing a 
power-law energy distribution and the Tsallis type anomalous thermodynamic rules. This model describes a non- 
equilibrium situation or a non-ergodic stationary state: some correlations are not suppressed exponentially, one 
particle distributions show a long tail in energy and hence in transverse momenta. 

This is exactly the case which has been observed in relativistic heavy ion collisions at RHIC at high transverse 
momenta. The Tsallis QGP is able to offer a temperature and energy per particle value which are near to the ones 
fitted by hadronic thermal models; a property not shared by traditional quark gluon plasma models. This way a 
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T (GeV) 

FIG. 5: Along the 1 GeV energy per particle line without bag constant in the T — /i plane the entropy per particle, S/N, is 
plotted. The constant power, c — 4.72732, line is solid, the constant energy parameter, b — 0.82728 GeV line is dashed-dotted, 
the traditional thermal line (c = cxd) is dashed. 




T (GeV) v^B (GeV) 

FIG. 6: Energy per particle in a Tsallis-QGP with adjusted bag constant to the p — line. Isochors (left) and isotherms (right) 
are shown. 

natural precursor state may be suspected at hadronization, as well as an explanation may be given for the value 
E/N Ri 1 GeV. Furthermore we have extended this estimate for finite baryochemical potential using Fermi-Tsallis 
and Bose-Tsallis statistics for quarks and gluons respectively. Also in this case one can come near to the hadronic 
thermal model curve with certain power fits, however, due to some indefiniteness in the treatment of anti-fermions 
the results are less convincing. 
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b = 1 .0325 GeV 
b = 2.0000 GeV 
b = infinity 
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FIG. 7: The E/N — 1 GeV line with p = adjusted bag constant on the T — /is plane {fiB = 3/i assumed). 

A possible solution to this problem is offered by introducing a monotonous re-scaling of Tsallis entropy, which leads 
to a straightforward correspondence between traditional and Tsallis distributions. The key to this is a logarithmic 
dispersion relation for the quasi-particles, which has an anomalous character, as we have discussed. This leads to an 
anomalous thermodynamics, in which now the anti-particles can be interpreted as holes in the particle distribution 
also for the extended, non-extensive Tsallis distribution. We have repeated the numerical analysis of the energy per 
particle and the spectra for this case, and realized that the basic qualitative features do not change. 

Our final conclusion from these studies is that in order to interpret experimental observations on the bulk particle 
spectra in relativistic heavy-ion collisions we need to use at last generalized thermodynamical concepts instead of 
the traditional ones. In particular only this approach allows an interpretation of a spectral slope, a high tail with 
power-law and the fitted E/N « 1 GeV energy per particle on the QGP side at the same time. 

It remains to clarify what mechanism leads to an anomalous statistics in these heavy ion collisions and to explain 
the very power. This is a task on its own right, here we just try to sketch some ideas about a possible explanation. 
Considering non-ergodic, but stationary solution of a Boltzmann equation for partons, special energy dependent cross 
sections may lead to an approximate validity of simple recombination rules |^ l25l l45|. The power-law distribution 
recombines in this case to distributions closer and closer to the familiar exponential [Ig. Analog to the "canonical 
suppression" , which considers effects due to a finite volume on the grand canonical distribution p^o . , here we aim 
to consider finite-time effects leading to partially-equilibrated distributions. 
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